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Abstract 

One-loop two-, three- and four-point scalar functions are analytically integrated di- 
rectly such that they are expressed in terms of Lauricella's hypergeometric function Fd- 
For two- and three-point functions, exact expressions are obtained with arbitrary com- 
bination of kinematic and mass parameters in arbitrary space-time dimension. Four- 
point function is expressed in terms of Fd up to the finite part in the expansion around 
f> |- 4-dimensional space-time with arbitrary combination of kinematic and mass parame- 

ters. Since the location of the possible singularities of Fd is known, information about 
the stabilities in the numerical calculation is obtained. We have developed a numerical 
library calculating Fd around 4-dimensional space-time. The numerical values for IR 
divergent cases of four-point functions in massless QCD are calculated and agreed with 
^ , golem95 package. 
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! 1 Introduction 

O 

One-loop calculation in perturbative field theory is a well-established and theoretically clear 
method. However, highly accurate numerical calculations of one-loop amplitudes are not a 
trivial problem. There appear many kinematic parameters including particle masses. It is 
k>( ■ not easy to keep numerical accuracy for every points in the multi-dimensional parameter 

; i \ space. Since the numerical behavior of a function is related to its analytic properties, it is 

important to pursue stable analytical expressions for numerical calculations. 

Feynman amplitudes are expected to be a kind of hypergeometric function II ]. One-loop 
integrals are explicitly expressed in terms of hypergeometric functions using several differ- 
ent methods: Mellin-Bares transformation, algebraic relations and power series expansion 
[UE1IUIS]. See also [5] for the references. In this article, we show that one-loop two-, three- 
and four-point functions are directly integrated with Gauss' function F, Appell's function 
Fx and Lauricella's i*£>[7]. Among these function, Fd includes other functions as special 
cases. Since the location of possible singularities of Fd is known[5|, we can point out dan- 
gerous combinations of parameters and to get information about where and how numerical 
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cancellation may occur. For two- or three-point functions, we show the integrals are ex- 
actly expressed in terms of Fd for any values of kinematic parameters in any space-time 
dimensions. For four-point case, we could not integrate exactly except some special cases. 
However, around 4-dimensional space-time, it is expressed in terms of Fd up to the finite 
order. 

We have developed a program package which calculate Fd for the necessary combinations 
of parameters for IR divergent case in massless QCD. We have compared our numerical result 
with golem95 package [TP]. 



2 Two-point point function 

Let's first consider two-point scalar function as a simple example of the usage of hypergeo- 
metric function. Two-point function is defined by: 



/>oo />oo 

4 a) = / dxi \ dx 2 S{1 - xx - x 2 ) T> a , 
Jo Jo 



(1) 

D = —p 2 xix 2 + m\x\ + m\x2 — ie, 

where, a is a number depending on the dimension of the space-time. It is noted that 
calculation with arbitrary mass parameters is useful even for massless cases, since tensor 
integrations are obtained by differentiating scaler one in terms of mass parameters. 
Integrating once with the 8 function, we obtain 



x \ x 
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(2) 



2p 2 ' v " 

This integral is nothing but a special case of Appell's Fi, whose integral form is: 

F 1 {a,b,b'yc ] y,z) = —^- j x a -\l - x)™-\\ - yx)-\l - zx)- b ' dx. (3) 



r(o)r(c-o) jo 

Thus we have 

4 a) = (^irm, ~a; 2; -ip, — ) - (m?)«fi(l, ~a, ~a; 2; — Lp -J— ). (4) 

The last equality is obtained by changing integration variable form x to y — 1 — x. This 
equality is considered as an identity of F±. 

Function F\ reduces to Gauss' hypergeometric function F for the case of b' = and is a 
special case of more general function, called Lauricella's Fd, whose integral form is: 



F D {aM,--- ,b n ;c; Zl ,--- , z n ) = ^ / x a -\l - z^" " 1 - z iX )- bi dx. 

I {a)l (c — a) 



■ / .T a - 1 (i-.T) c - a - 1 TT(i 

J ° i= i 



(5) 



This function Fd will be used for vertex and box integration. The location of possible 
singularities of Fd is limited to z% = 0,1, 00 and Zj — Zfc (j 7^ fc) and Fd is smooth except 
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these special cases [5]. The power series expansion and differential equation of this function 
is known. Many of identities of F can be generalized to Fr> . We especially use the following 
identity, which can easily be confirmed from Eq. ([5]) : 



z p-i {1 _ z y-i JJ (1 _ x . z yU = ^-F D (p, (b t ), 1 - q;p+ 1; ( Xi z),z) . (6) 

This identity implies that any product of linear factors with arbitrary power is integrated 
by F D . 

Now going back to Eq.(j4j). Two-point function may be singular when 7 ± = 0, 1, oo 
or 7 + = 7~. These cases correspond to massless particles, p 2 — 0, and on the threshold. 
Let us examine the case of rn| = and m 2 ,p 2 ^ for an example, where 7 _ = and 
7 + = (p 2 — mf ) jp 2 . Although Eq.@ is not well-defined for the limit of m-i — > 0, we can use 
alternative representation obtained with another identity of F±: 

= -l—( m 2 ) a F{a + l,-a- 1 a + 2;^- ^ 



a + l 7 - 7+ 

+ ^—^-{ml) a F{a + 1, -a; a + 2; * " 7 _ ). 

This representation is regular for the limit of rri2 — > under the condition of 5ia > 0. Then 
expanding around 4-dimensional space-time ( a = — e — > +0 ), F\ is reduced to F and 
then to log and Li2. Two limiting processes m,2 — > and a — > +0 do not commute. This 
property corresponds to the fact that the analytic result in 4-dimension does not reduces to 
the massless one by taking the simple limit ttt-2 — > 0. 

This example shows that Eq. @ can be used as a unified representation for both massive 
and massless cases. When values of mass parameters are specified, we select a suitable 
representation with appropriate limit. We want to select it not on a notebook but in a 
numerical library at the time of numerical calculations. With this library, the main program 
will be general to cover various cases of parameters. The problem of numerical instability 
will be confined into the numerical calculation method of Fd ■ 



3 Three-point function 

Three-point function is defined by: 



I L 3 a) = / dx l dx 2 V a . (8) 

J X 1 ,X2>0,X 1 +X2<1 

where, T> is a quadratic form of x\ and xi- We apply the projective transformation as shown 
by Ref. [9 . The quadratic term of Xi is eliminated by changing variables (xi, x-z) — ¥ (x2, z = 
xi + rx2) with adjusting of the value of r. Since T> is now linear in X2, integration is trivial 
for Xi- The resulting integration becomes the form: 



J 3 Q) « / ^TT-T dz - ( 9 ) 



az + b 

As T> is expressed as a product of linear factors of z, we obtain: 
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Using Eq.©, this is immediately integrated by The integration domain becomes slightly 
complicated after the projective transformation. It is handled systematically with exterior 
derivative and Stokes' theorem. The result takes the following form: 

^3 =— Th- - — -r^F D (l,l,-a-l,-a-l;2;- — ,— ,— ) (11) 

where, is the value of I? at the corner of the integration domain and d k j is brought from 
the parameterization of the boundary after the projective transformation. For the limit of 
4-dimcnsional space-time, the above expression reduces to the usual analytic representation 
as Fd reduces to Fjj — > F\ — » F — > log and poly-logarithmic functions. 




4 Four-point function 

Four-point function is written by: 

I<f ] =l 2>«. (12) 

Ho V J 

where I? is a homogeneous quadratic form of Xj. After using 5 function, there left three 
integration variables. 

We apply projective transformations twice. After the first transformation, we integrate 
once using the following identity: 

1 T> a+1 

V a dy = — — , (13) 

,y a + 1 d v V ' v ; 

where T> is quadratic in terms of remaining two variables, while d y T> is linear. 

The projective transformation is applied once more. T> becomes a linear function of a 
new variable z. It is possible to select the variable z by shifting and rescaling such that 
z oc T> and 1 — z cx d y T>. Integral is calculated with the following formula (special case of 
Eq. ©): 

1 — z dz 

where F is Gauss' hypergeometric function. After the second integration we obtain: 

4 Q) =E E Cfef/ [ 9k + h k (e k )]dy ke , (15) 



fe=l e=i,e=£k 



3 fc +1 f d k V k \ a+2 d k T> k . 



(a + l)(« + 2 K «V « J a + 2,a + 3; — ), (16) 

(k) 

where D k and e k are quadratic form of integration variable y k i , d k and Q are brought by 
projective transformation, and L k i is a line segment of the last integration. Function h k of 
e k is arbitrary and is produced as an integration constant. 

In order to handle F in the integrand, we use partial integration method. Using recursion 
relation of F, g k is expressed by: 

1 e. k / d k V k 



9k (a + l)(a + 2)d«+ 2 V e k 



(17) 

1 e k ^ f d k V k \ d k V k 

F(l,a + 3,a + 4; ) 



(a + l)(a + 3)d«+ 2 \ e k J v ' ' ' e k 
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Factor e? +1 is integrated by the following relations: 



e k {x) = e k (w 5 - x)(w 6 - x), 



„a+l 



(X) 



df(x) 
dx 1 



i 



(,;) 



a + 2 e k (w 6 - w 5 ) 
After the partial integration, we obtain: 

4 4 



1 - 2F(-a - 2, a + 2, a + 3; — — -) 

wq — x 



r(«) _ 
1 4 — 

Jl = 
J 2 = 
■h = 



k=l,k=£m £=l,i=ik 



1 



1 v k a+2 



(a + l)(a + 2) 7o e fc 
1 



1 



(a + l)(a + 3)d«+ 2 



d k V k 



Q + 3 



F(l,a + 3,a + 4; 



e k 



x=0 



1 1 



1 rl a + 2 \ 

1 "fc Jo 



d k V k 
e k 



a+2 



d 1 

y log 

dx 



e k - d k V k 



( : 'k 



dx. 



Here we have used the following identity: 



4-i?(x) a+3 F(l, a + 3, a + 4; = -(a + 3)i?(x) a+2 4- logfl - R(x)]. 

dx dx 



(18) 



(19) 



(20) 



Ji is integrable with Fjj, since the integrand is expressed by a product of power of linear 
factors. Ji is a product of F. The problem left is function / in J3. 

Investigating the limit to the 4-dimensional space-time a = — 2 — e — > —2, one will 
confirm that integration in J3 does not produce new poles of 1/e. So we can expand in 
terms of e in the integrand. 



F(e,-e,l-e;z) = l + 0(e 2 ). 



■h 



e(l + e) e fc (u; 6 - w 5 ) J 



V k a+2 ^\02 



d 

dx 



c-k — d k V k 



('k 



dx + 0(e) 



(21) 
(22) 



Since factor dlog j dx is expressed by a sum of inverse of linear term of x, J3 is expressed by 
Fd up to the finite order. 



5 Sample numerical calculation 

Since F, F\ and Fp have many parameters and variables, it is hard to construct numerical 
package to calculate for all cases. However, we need numerical values only for some special 
combination of parameters for our purpose. 

We have tried sample numerical calculations of one-loop box tensor integration for the 
cases of Ref. g]. That is: 

• All particles are massless. 

• At least one external particle is on-shell (p\ = 0). 

• Calculate up to O(e ). 
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• Calculate tensor integrations up to rank = 4. 



There appear IR divergences, which are represented by the poles of 1/e. 

For these cases, integration becomes simpler by using variable transformation as de- 
scribed in [3] . We have obtained an exact analytic representation with Fd for the cases of 4 
or 3 on-shell particles and "easy case" of 2 on-shell particles (diagonal external particles of 
the box diagram are on-shell). However, for "hard case" of 2 on-shell particles (two adjacent 
external particles of the box diagram are on-shell) and 3 on-shell case, we need to expand 
in terms of e. 

It is noted that representations of n-point functions are not necessarily numerically stable 
when they are written with poly-logarithmic functions. For example, there appears F(l, to — 
e, m + 1 — e; z) (to > 1) in the tensor integrations, which is regular around z ~ 0. When this 
function is expanded in terms of e using identities of F, the following combination of terms 
appears: 



When Lifc + i is expanded around z ~ 0, the first m — 1 terms of the power series cancel 
out with the second term in the brackets. Factor z m is factored out from the resulting 
terms in the bracket, and cancels with the denominator. If these terms are scattered into 
a long expression, it is not an easy problem to control these numerical cancellations and 
the singular behavior of the denominator. However, when we keep the original form of 
F(l, to — e, to + 1 — e; z), the problem is immediately solved; a simple power series calculation 
of F around z = produces stable result. This problem is caused by the expansion made 
in order to express these functions in terms of Life. 

We have developed a numerical library of hypergeometric function Fd for necessary com- 
bination of parameters for our sample calculations. Our library is designed in the following 
way: 

1. Entry points of subroutines are Fd or F. 

2. These subroutines return an array of coefficients of 1/e 2 , 1/e, 1, e, • • • up to necessary 



3. Inside of the subroutines, appropriate identities or calculation methods are selected in 
looking at the values of parameters and variables. 

Two programs are prepared for the numerical calculations of tensor integrations: 

• "program-1": calculation with numerical library of F, F\ and Fd- 

• "program-2" : calculation with the following numerical integration: 

— The first two integrations are calculated analytically. 

— Coefficients of 1/e 2 , 1/e 1 , l/e° are extracted and expressed by one-dimensional 
integrations. 

— The last integration is calculated numerically (Romberg method). 




k 



m—l 



order. 
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We have compared the numerical results among program-1, program-2 and golem95 
package up to rank = 4 at 7560 different values of the parameters: 

p\ = Q, p| = Q,±50, p^ = 0,±55, pl = 0,±60, 
s = ±200, t = ±123, 

Ui = 0, 1, 2, 3, 4, m < 4 (rank of tensor integration) 

i 

The results of the maximal differences among methods are shown by Table [TJ It shows that 
the accuracy of the library seems similar to golem95 package. 



calculation method maximal difference 



program- 1(d) 


program-2 (d) 


7.65 x 10" 


-7 


program- 1(d) 


golem95 (d) 


9.13 x 10" 


-10 


program- 1(d) 


program-1 (q) 


3.98 x 10" 


-10 


golem95(d) 


golem95(q) 


5.17 x 10- 


-10 


program-1 (q) 


golem95 (q) 


1.38 x 10" 


-18 



Table 1: Maximal differences among the calculation method. Differences are measured by 
the distance on the complex plane, (d) and (q) stand for double and quadruple recision 
respectively. 



6 Summary 

Two- and three-point functions are expressed in terms of Fd, exactly for any combination 
of physical parameters in any space-time dimensions. Four-point functions are expressed 
with Fd, up to 0(e ) for any combination of physical parameters. A program library of 
F and Fd is developed applicable for sample numerical calculations for massless QCD box 
with IR divergences. The results agree with golem95 package. Four-point function seems 
not to be integrated with Fd ■ In order to express general four-point function, more general 
hypergeometric functions will be needed as described in [31 [3] and [TT] for more general 
cases. 
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